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ABSTRACT 

We explore whether the rest-frame near-ultraviolet spectral region, observable in 
high-redshift galaxies via optical spectroscopy, contains sufficient information to allow 
the degeneracy between age and metallicity to be lifted. We do this by first testing 
the ability of evolutionary synthesis models to reclaim the correct metallicity when 
fitted to the near-ultraviolet spectra of F stars of known (sub-solar and super-solar) 
metallicity. F stars are of particular interest because the rest-frame near-ultraviolet 
spectra of the oldest known elliptical galaxies at z > 1 appear to be dominated by F 
stars near to the main-sequence turnoff. 

We find that, in the case of the F stars, where the HST ultraviolet spectra have 
high signahnoise, model-fitting with metallicity allowed to vary as a free parameter is 
rather successful at deriving the correct metallicity. As a result, the estimated turnoff 
ages of these stars yielded by the model fitting are well constrained. Encouraged by 
this we have fitted these same variable-metallicity models to the deep, optical spectra 
of the z ~ 1.5 mJy radio galaxies 53W091 and 53W069 obtained with the Keck tele- 
scope. While the age and metallicity are not so easily constrained for these galaxies, 
we find that even when metallicity is allowed as a free parameter, the best estimates 
of their ages are still > 3 Gyr, with ages younger than 2 Gyr now strongly excluded. 
Furthermore, we find that a search of the entire parameter space of metallicity and 
star formation history using MOPED leads to the same conclusion. Our results there- 
fore continue to argue strongly against an Einstein-de Sitter universe, and favour a 
A-dominated universe in which star formation in at least these particular elliptical 
galaxies was completed somewhere in the redshift range z = 3 — 5. 

Key words: galaxies: evolution -- galaxies: stellar content - ultraviolet: galaxies - 
ultraviolet: stars 



1 INTRODUCTION 

Deep optical spectroscopy with 8 — 10-metre class telescopes 
can overcome many of the problems associated with estimat- 
ing the ages of high-redshift galaxies from cruder broad-band 
data. In particular, as demonstrated by Dunlop et al. (1996), 
Spinrad et al. (1997) and Dunlop (1999), high-quality spec- 
tral data not only provide reassurance that one is observing 
starlight, but allow reddening-independent age estimation 
based on the strength of specific spectral features. These 
observational advances have greatly enhanced the prospects 
for reliable age-dating of high-redshift stellar populations. 
Consequently, recent years have seen renewed interest in 
clarifying and attempting to resolve the remaining sources 
of uncertainty. As is often the case when an area of high- 
redshift astronomy finally enters the regime of quantitative 
astrophysics, the dominant sources of uncertainty transpire 



to be problems that have dogged work at lower redshift for 
many years. In this case the main outstanding issues are 
disagreements over the details of stellar evolution, and the 
well-known degeneracy between derived age and metallicity 
(Worthey 1994). 

Of course the potential severity of both, and especially 
the former of these effects is expected to be a function of 
age and observed spectral region. In particular, as described 
by Magris & Bruzual (1993), and reaffirmed in the study of 
the red, z ~ 1.5 galaxy LBDS 53W091 (Dunlop et al. 1996; 
Spinrad et al. 1997), it is in the rest-frame near-ultraviolet 
(near-UV), at ages less than 5 Gyr, that disagreements be- 
tween the predictions of alternative models of spectral evo- 
lution should be minimised. This is because for such ages 
it is expected that the near-UV light should be dominated 
by the output of stars near the turnoff point of the 'well- 
understood' main sequence (e.g. Nolan et al. 2001). How- 
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ever, the ensuing controversy over the precise age of 53W091 
(Bruzual & Magris 1998; Yi et aL 2000), and the resulting 
confusion over the wider implications for cosmology in gen- 
eral, have highhghted the need for a closer inspection of the 
uncertainties which continue to afflict galaxy age-estimation. 

In an attempt to clarify the key issues we have now 
begun to study both the extent to which various evolution- 
ary synthesis models differ in their evolutionary predictions, 
and the extent to which rest-frame near-UV data may allow 
the age-metallicity degeneracy to be lifted. In Nolan et al. 
(2001), we checked the time evolution of the solax-metallicity 
(Zq) population synthesis models of several authors, by 
comparing the ability of the main-soquoncc components of 
these models to reproduce the near-UV solar spectrum of 
the Sun at its anticipated main-sequence turn-off age. 

In this paper we explore whether, given data of suf- 
ficient quality, variable-metallicity stellar population mod- 
els can extract the correct mctallicity of a population sim- 
ply from fitting to the near-UV spectral region accessible 
through optical spectroscopy of high-redshift galaxies. 

We test whether the UV spectral region contains suffi- 
cient information to allow reliable recovery of the metallicity 
in two stages. Firstly, the single stellar atmosphere model 
spectra of Jimenez et al. (1998) are fitted to the HST ul- 
traviolet spectra of two well-studied, non-solar metallicity, 
main-sequence stars, at a resolution A/AA = 30,000 (Heap 
et al., 1998a), and a typical signal-to-noise per pixel of 50, 
before re-binning (Heap et al., 1998b). These are kindly sup- 
plied by Sally Heap (Heap, private communication). Sec- 
ondly, variable- and mixed-metallicity versions of the stellar 
population models of Jimenez ot al.(1998) arc fitted to the 
same two stellar spectra. These stars are HR 4683, an F4V 
star, with sub-solar metallicity, and HR 4688, an F8V star 
with super-solar metallicity (Edvardsson et al., 1993). The 
motivation for using the spectra of these particular stars is 
that the rest-frame near-UV spectra of the oldest known 
galaxies at 2; ~ 1.5 (e.g. 53W091 and 53W069) mimic the 
spectra of main-sequence F stars (F5V - F9V; Dunlop 1999). 

This is, of course, a rather one-sided test, but it is still 
potentially very instructive. Specifically, if the single stel- 
lar atmosphere models cannot extract the correct mctallic- 
ity from the high-quality near-UV data available for these 
F-stars, then it is probably safe to assume that lifting the 
age-metallicity degeneracy in high-redshift galaxies is cur- 
rently impossible. Conversely, if the single atmosphere mod- 
els can yield the correct metallicity for these F-stars, then 
that provides some hope that the age-metallicity degener- 
acy can be lifted in high-redshift galaxies, given near-UV 
data of suflicicnt quality. In addition, if the variable- and 
mixed-metallicity stellar population model results are con- 
sistent with the results from fitting the single stellax atmo- 
spheres, then we may be confident that, at young ages, main- 
sequence stars do indeed dominate the UV flux of stellar 
populations, that the time-evolution of the non-solar metal- 
licity population models is correct, and that we can con- 
strain metallicity even with the greater flexibility allowed 
in the mixed-metallicity population models. We may then 
be confldent that ages and metallicities recovered for the 
high-redshift galaxies arc reliable. 

However, such optimism must be tempered by the 
knowledge that the spectrum of a high-redshift galaxy is un- 
likely to be as clean as that of the test-case F-stars, given the 



obvious likelihood of a mixture of metallicities and ages in a 
substantial stellar population. Nevertheless, as we report in 
this paper, the results of fitting the variable/mixed metallic- 
ity models of Jimenez et al. (1998) to the near-UV spectra of 
the two F-stars are sufficiently encouraging to motivate the 
re-fitting of the Keck spectra of LBDS 53W091 and 53W069 
with metallicity allowed to vary as a free parameter. 

The layout of the paper is as follows. In section 2 we de- 
scribe the models we have used, and the procedures by which 
they have been fitted to the data. In section 3 we present 
the results of fitting single, variable metallicity, model stellar 
atmosphere spectra, single (but variable) metallicity stellar 
population models, and mixed-metallicity stellar population 
models (comprising a free mix of seven different metallici- 
ties) to the HST near-UV spectra of the F stars. Then, in 
section 4, we present the results of fitting the same stel- 
lar population models to the optical (rest-frame near-UV) 
spectra of the z ~ 1.5 galaxies LBDS 53W091 and 53W069, 
and discuss the extent to which such modelling allows the 
age-metallicity degeneracy to be lifted. Finally, as a check 
on whether our assumption of single burst models is appro- 
priate, we have also fitted the spectra of these two galaxies 
using MOPED (Heavens et al., 2000), which allows the star 
formation history to be unconstrained. Our conclusions, in- 
cluding revised constraints on the ages of these red galaxies, 
are summarised in section 5. 



2 MODELS AND MODEL FITTING 
2.1 The models 

The stellar population evolutionary synthesis models 
utilised in this paper are the instantaneous starburst models 
of Jimenez et al. (1998), and the model single stellar atmo- 
spheres are those used in creating these population mod- 
els. The reliability and consistency of the solar metallicity 
(Z0= 0.02) version of these population models were con- 
firmed by Nolan et al. (2001). We now also possess non-solar 
metallicity models, for both the stellar population models 
and the single stellar atmosphere models, with metallicities 
0.01 Zq, 0.2 Zq, 0.5 Zq, 1.5 Zq, 2.5Zoand 5.0 Zq. 

The stellar interior library of the Jimenez et al. (1998) 
models is computed using the latest version of the JM- 
STAR15 code (Jimenez et al. 1996, Jimenez & McDonald 
1996). This uses the OPAL95 opacities (Iglcsias & Rogers 
1996) for stellar temperatures ~ 6000K, and Alexander's 
opacities for temperatures < 6000K. The mixing- length the- 
ory of Baker & Temensvary (1966), modified to allow for the 
limit of optically thin convective elements (Mihalas 1978), 
is used for the treatment of convective energy transport. 
The model grid is calculated for a range of metallicities, Z 
= 0.0002 - 0.1, with solar metallicity, Z© = 0.02, and the 
solar helium fraction, Y© = 0.28. The mixing-length param- 
eter, a, is calibrated by a comparison of the computed solar 
model with the observed effective temperature, Teff, age, 
luminosity and radius of the Sun. In addition, a consistency 
check was carried out, using the position of the red giant 
branch (RGB) in globular clusters with known metallicity 
(Jimenez et al. 1996). 

Stellar masses range from 0.1 to 120 M©, with step 
size, A M = 0.05 M©, from 0.1 to 3 M©, and A M = 1 M© 



© 0000 RAS, MNRAS 000, 000-000 



F stars, metallicity, and ages of red galaxies 3 



thereafter. The models used in this work have a Salpctcr 
IMF (Salpeter 1955). The tracks are evolved all the way 
from the contracting Hayashi track to the planetary nebula 
(PN) phase, except for the most massive stars, where the 
evolution is followed until carbon ignition. All tracks are 
computed using solar-scaled abundances. 

Stellar evolutionary tracks axe computed using the 
boundary conditions set by the theoretical photospheric 
models calculated by Jimenez et al. (1998). Therefore, for 
each calculated point on the track, the associated stellar 
spectrum is known, and there is no need to introduce the 
uncertainties which arise in calibrating real stellar spectra 
to theoretical isochrones. Semi-convection arises as a natural 
consequence of this modelling method. 

JMSTAR15 computes the whole evolution of a star in 
a single run, including the helium core flash, through the 
thermally-pulsing asymptotic giant branch (TPAGB), eject- 
ing a PN at the end. This is possible due to the adaptive 
mesh of grid points, which allows very rapid stages of evo- 
lution to be modelled in detail. The SSP code includes a 
semi-empirical algorithm for computing the evolution of the 
RGB, horizontal branch (HE) and asymptotic giant branch 
(AGE). This is more realistic than a purely numerical code, 
and faster, but allows the same degree of accuracy. The AGE 
was modelled as in j0rgensen (1991). 

A semi-empirical algorithm is also used in the calcu- 
lation of mass-loss. The fast computation of post-main se- 
quence (post-MS) stages allowed several values of the mass- 
loss efficiency parameter, {r} in Reimers' formula, Reimers 
1975), and mixing- length parameter, a, to be investigated. 
A realistic mean value of rj is then found by matching the 
observed HE mass distribution (Jimenez et al. 1995). 

The stellar spectra arc those of the Kurucz (1992) li- 
brary of theoretical stellar atmosphere models for stars with 
Te// > 6000 K. For cooler stars, theoretical photospheres 
are computed by Jimenez ct al. (1998), using the improved 
(Helling et al. 1996) MARCS code (Gustafsson ct al. 1975), 
which includes molecular opacities. 

2.2 The fitting processes 

Firstly, we have explored the results of simply fitting each 
of the seven single stellar atmosphere models in turn to the 
F star spectra, at all possible ages, to identify the best- 
fitting single value of metallicity, and the resulting inferred 
age. Secondly, the same process is carried out for the seven 
single-metallicity stellar population models, for both the F 
star spectra and the high-redshift galaxy spectra. Finally, 
we have combined the seven single-metallicity stellar pop- 
ulation models to produce a composite model (of single 
age) in which mean metallicity is allowed to vary by vary- 
ing the fractional contributions made by each of the seven 
different metallicity components. We have fitted the mixed- 
metallicity model to the spectra of the F stars, as well as the 
high redshift galaxies' spectra, to confirm that fitting this 
model allows us to reclaim the correct metallicity. 

2.3 Fitting the single stellar atmosphere and 
single-metalhcity population models 

The best fit was determined for each single-metallicity model 
(single-atmosphere and stellar population), i.e 0.01 Z0, 



0.2 Z0, 0.5 Zq, Z0, 1.5 Z0, 2.5Z0and 5.0 Z0, in the follow- 
ing way. First, the data were re-binned to the same spectral 
resolution as the models, with the new flux calculated as 
the mean flux per unit wavelength in the new bin. Both 
the rc-binncd and model fluxes were then normalised to a 
mean flux per unit wavelength of unity across the wave- 
length range 3000 — 3500A for the galaxy fits, and across 
the range 3000 — 3117A for the F-star fits (because of the 
shorter wavelength range of the F-star spectra). The model 
age and normalisation were then varied as free parameters 
until X was minimised. The wavelength range 2445-2465 A 
was excluded from the fits to all objects. When included, this 
feature contributed to 1/3 of the total for both the F 
stars. As this was felt to unfairly bias the fitting process 
for these objects, the area was masked from the fit for all 
the objects, and may indicate a failure of the models in this 
range. 

In order to produce the contour plots for the single- 
metallicity stellar atmosphere and stellar population model 
fits in sections 3 and 4, the results were refined by carrying 
out two one-dimensional cubic spline interpolations, one for 
age, and one for metallicity, between the calculated points 
on the age-metallicity grid. 



2.4 Fitting mixed-metallicity models 

A mixed-metallicity model was constructed from the 
0.01 Z0, 0.2 Z0, 0.5 Z0, Z0, 1.5 Z0 2.5 Z0 and 5.0 Z0 stellar 
population models of Jimenez et al. (1998). The relative con- 
tributions of the different metallicity models were allowed as 
free parameters. In this section, we assume that all the com- 
ponent populations were formed at the same time, so that 
the one other free parameter in the fit was the age of the 
complete population. In section 4.3, we lift this constraint 
and explore what happens when the stellar populations are 
formed at arbitrary times. 

The mixed-metallicity model flux was built from the 
unnormalised single-metallicity model fluxes, and then nor- 
malised to a mean flux per unit wavelength of unity in the 
same regions as for the single-metallicity model fits, so that 

7 

F7z,x,t = consty^^Xzifzi,x,t 

i=l 

where F-jz,\,t is the mixed, seven-metallicity model flux per 

unit wavelength in the bin centered on wavelength A at age, 
t Gyr; fzi,\,t is the flux per unit wavelength A at age, t Gyr, 
for the single metallicity model with Z = Zi, where i = 1,7, 
and Xzi is the fractional contribution (by baryonic mass) of 
the model with Z = Zi to the total population. 

The best-fitting age, normalisation and the value of each 
Xz^ were determined by minimisation, with the wave- 
length range 2445 — 2465 A masked from the fit, as it was in 
the single stellar atmosphere and single stellar population 
fits. The entire parameter space of the seven-dimensional 
age-metallicity hyper-cube was searched, so that the best 
fits quoted in sections 3 and 4 are those parameter values 
which correspond to the point on the age-metallicity grid 
with the minimum x^- 
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Figure 1. Top: the best-fitting, single stellar atmosphere model 
(thin line, 0.2 Zq model at 4 Gyr) to the spectrum of the F4V 
star, HR 4683 (thick line). Bottom: residuals. Where {f\data < 
ixmodel), the quantity plotted is -{{{xdata - hmodel) / '^)'^ ■ The 
central wavelengths of the mid-UV spectral features identified 
by Fanelli ct al. (1992), lying within this wavelength range, are 
marked on the spectra with dashed lines. See §3.1 for discussion, 
and table 2 for values of ■ 
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Figure 2. Top: the best-fitting, single stellar atmosphere model 
(thin line, 2.5 model at 1.5 Gyr) to the spectrum of the F8V 
star, HR 4688 (thick line). Bottom: x^ residuals. Where {i\data < 
hmodel), the quantity plotted is -{{hdata - ^Xmodel) I <^)'^ ■ The 
central wavelengths of the mid-UV spectral features identified 
by Fanelli et al. (1992), lying within this wavelength range, are 
marked on the spectra with dashed lines. See §3.1 for discussion, 
and table 1 for values of x?- 



2.5 Errors 

The errors on the binned spectral data points for the LBDS 
gala^xies have been derived from the propagation of the orig- 
inal observational errors in the Keck optical spectroscopy 
(Dey et al., 2002, in preparation). Using these errors, statis- 
tically acceptable fits were achieved. For the F-star spectra, 
the error adopted for each rc-binncd data point was simply 
the standard error in the mean for that bin. For HR 4683 
(the low-metallicity star), a formally acceptable fit was ob- 
tained using these errors. However, for HR 4688 (the high- 
metallicity star), this did not prove possible. There are a 
number of reasons why this might occur (e.g. inadequacy of 
the super-solar models — e.g. Kotoneva, Flynn & Jimenez 
(1998). 



2.6 Marginalising the mixed-metallicity results 

The contours of constant relative likelihood discussed in §3.3 
and 4.2 arc plotted with respect to mean motallicity and age. 
As it is possible for more than one combination of single- 
metallicity models to result in the same mean metallicity, the 
metallicity distribution was marginalised before plotting. 



3 COMPARISON WITH F-STAR SPECTRA 
3.1 Single stellar atmospheres 

In table 1 and figures 1 and 2, the results of fitting the seven 
single-atmosphere models, with metallicities 0.01 Z0, 0.2 Zq, 
0.5 Zq, Zq, 1.5 Zq, 2.5 Zoand 5.0 Zq to the HST spectra of 
the two F stars are presented. The best-fit single-star model 
to the F4V star spectrum has a metallicity of 0.2 Z©, and 
an age of 4 Gyr. For the F8V star, the best-fit single-star 
model has an age of 1.5 Gyr, with a metallicity 2.5 Zq. 

Edvardsson ct al. (1993) estimated the relative iron 
abundances ([Fe/H] = Z/Zq, where Zq = 0.02) of these 
two stars by fitting model stellar atmospheres to spectral 
lines observed at optical wavelengths. They estimated the 
uncertainty in the derived metallicities to be, at most, 0.10 
dex. Based on these optical lines, they found the F4V star 
to have a metallicity of 0.29 Zq, and the F8V star to have a 
metallicity of 1.62 Zq. 

Edvardsson and coworkers also derived ages for these 
stars from fits in the Te// — log g plane to VandcnBerg (1985) 
isochroncs. They found an age of 4.4 Gyr for the F4 dwarf, 
HR 4683, and an age of 2.1 Gyr for the F8 dwarf, HR 4688. 
However, they expect the uncertainties in the relative ages 
of the stars in their sample to bo ~ 25%, with the possibility 
that the absolute errors are larger. 

Our results from fitting the single stellar atmosphere 
models are in excellent agreement with these, entirely inde- 
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Figure 3. Top: the best-fitting, single-metaUicity stellar popula- 
tion model (thin line, 0.2 Zq model at 3 Gyr) to the spectrum 
of the F4V star, HR 4683 (thick line). Bottom: contours of con- 
stant reduced (x5) ^ function of metallicity and age after 
interpolation. The contour levels contain 68%, 90% and 95.4% 
relative likelihood. The shaded regions represent the estimates of 
Edvardsson et al., with estimated errors ± 30% in age and ± 0.10 
dex in metallicity. residuals are shown in figure 5. See §3.2 for 
discussion, and table 2 for values of Xv- 
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Figure 4. Top: the best-fitting, single-metaUicity stellar popula- 
tion model (thin line, 1.5 Zq model at 3 Gyr) to the spectrum of 
the F8V star, HR 4688 (thick line). Bottom: contours of constant 

as a function of metallicity and age after interpolation. The 
contour levels contain 68%, 90% and 95.4% relative likelihood. 
The shaded regions represent the estimates of Edvardsson et al., 
with estimated errors ± 30% in age and ± 0.10 dex in metallicity. 

residuals are shown in figure 6. See §3.2 for discussion, and 
table 2 for values of xt ■ 



pendently obtained, age and metallicity estimates. Encour- 
aged that the UV spectral region contains sufficient data to 
determine both age and metallicity, we proceed to fit the 
stellar population models. 

3.2 Single- metallicity stellar population models 

Attempting to fit complete stellar population models to the 
spectra of single stars might seem like a strange thing to 
do. However, having established that the individual input 
stellar spectra seem to work well in the near-UV regime, 
extending the model fitting to complete stellar populations 
ofTers a means to test the extent to which the near-UV light 
from a complete stellar population is indeed dominated by 
the dwarfs near the MS turnofT. It is also important to test 
whether, with the addition of stars spanning a wide range in 
mass and evolutionary stage, the theoretical UV spectrum 
can still identify the correct metallicity of our adopted turn- 
off population (as represented by the single F-star spectrum) 
given data of this quality. Finally, expanding the parameter 
space still further to include a mix of stellar populations of 
different metallicity, it is of interest to test whether these 
stellar spectra still contain enough information to allow the 
age and metallicity to both be well constrained. 

In figures 3 and 4, we show the results of fitting each 
of the seven single-metaUicity stellar population models in 



turn to the HST spectra of the F stars HR 4683 and HR 
4688 respectively. The first plot in each figure shows the 
best-fitting model spectrum superimposed on the observed 
near-UV stellar spectrum. The second plot shows contours 
of constant reduced ixi) ^ function of metallicity and 
age. The contour levels contain 68%, 90% and 95.4% relative 
likelihood. The results are presented in table 2. 

The fitting process results in a minimum xt correspond- 
ing to an age of 3 Gyr, with a metallicity of 0.2 Z© for the 
F4V star, HR 4683, and an age of 3 Gyr, with a metallicity 
of 1.5 Zgfor the F8V star, HR 4688. 

Again, the metallicities yielded by our completely inde- 
pendent, near-UV model-fitting shown in figures 3 and 4 are 
in very good agreement with those previously derived from 
optical indices by Edvardsson et al. (1993); we have derived 
metallicities from the near-UV data which fall within the 
uncertainties in the optically-derived values (see figures 3 
and 4). 

It should be noted that the ages determined by Ed- 
vardsson and co-workers are the present ages of the stars, 
whereas fitting the stellar population models to the spectra 
should yield ages which correspond to the main-sequence 
lifetimes of the stars. For this reason it might be expected 
that our model-fitting should yield ages greater than those 
derived by Edvardsson and co-workers. However, the mod- 
els include integration down the main- sequence, as well as 
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object Z/Zq Best fit age/Gyr y^i, object Z/Zq Best fit age/Gyr 



F4V 



0.01 
0.2 
0.5 

1.0 
1.5 
2.5 
5.0 



4 

3 
2 
1.5 
1 

0.5 



2.47 
1.59 
1.73 

2.43 
2.71 
4.00 
8.31 



F4V 



0.01 
0.2 
0.5 
1.0 
1.5 
2.5 
5.0 



13 
3 
3 
1 
1 

0.7 
0.3 



10.5 
1.68 
1.72 

2.79 
4.48 
5.29 
8.45 



F8V 



0.01 
0.2 
0.5 
1.0 
1.5 
2.5 
5.0 



not on MS 
5 
3 
2.5 
2 
1.5 
1.5 



3.19 
3.21 
2.97 
2.58 
2.40 
3.20 



F8V 



0.01 
0.2 
0.5 
1.0 
1.5 
2.5 
5.0 



13 
13 
13 
4 
3 
1 

0.7 



21.7 
2.86 
2.55 
2.52 
2.35 
2.47 
3.76 



Table 1. The results of fitting the ncar-ultraviolct spectra of 
the F-stars using the suite of seven single-metallicity single- 
atmosphere stellar models. The best fits for each of the single- 
metallicity stellar atmosphere models corresponding to the mini- 
mum xt shown. 

post main-sequence contributions, both of which are likely to 
lead to an undcr-cstimatc of the true main-sequence turnoff 
age when fitted to the spectrum of a single star. Therefore, 
more important than the exact values of the derived ages 
for single stars is the fact that the age appears to be well 
constrained despite the availability of seven different metal- 
licities at each age. In addition, the results of these stellar 
population model fits are in good agreement with the results 
of §3.1, as one would expect if stellar populations axe indeed 
dominated by MSTO stars at ages < 5 Gyr. 

Those results therefore demonstrate that by fitting the 
single-metallicity models to the near-UV spectra of the F 
stars, we have apparently managed to simultaneously deter- 
mine both age and metallicity, recovering well-constrained 
and plausible metallicities for these two stars. 

3.2.1 Age / metallicity determinators 

This success leads naturally to the question of whether any 
specific spectral features in this near-UV regime can be iden- 
tified as primarily responsible for allow the metallicity and 
age to be so efi^ectively determined in this analysis. To ex- 
plore this we have derived residual plots for a range of 
single-metallicity stellar population models fitted to the ob- 
served stellar spectra, which are shown in figures 5 and 6. 
The central plot is the residual from the best over-all fit. 
The metallicity of the model fitted increases downwards, and 
the age increases from left to right. Where the (normalised) 
observed flux is less than the (normalised) model flux, the 
quantity plotted is -{hmodei - hdata/cfT- The feature at 
2450 A has been masked out in all fits. 

In fact, from these figures it is hard to identify any par- 



LBDS 53W069 0.01 13 2.69 

0.2 13 0.92 

0.5 9 0.92 

1.0 4 1.04 

1.5 3 1.02 

2.5 1 1.10 

5.0 0.7 1.26 



LBDS 53W091 0.01 11 3.17 

0.2 10 1.26 

0.5 8 1.32 

1.0 3 1.21 

1.5 3 1.30 

2.5 1 1.37 

5.0 0.7 1.51 



Table 2. The results of fitting the near-ultraviolet spectra of 
the F-stars and the LBDS galaxies using the suite of seven single- 
metallicity stellar population models. The best fits for eaich. of the 
single-metallicity models are given. The best fitting models in the 

case of the F4V star, 53W069 and 53W091 arc acceptable fits - 
very acceptable in the case of the two galaxies. However, the best 
fitting single metallicity model to the F8V star is not formally 
acceptable without further relaxation of the assumed errors. See 
§3.2 for discussion. 



ticular absorption line as being primarily an indicator of age 

or of metallicity. However, there is a clear suggestion that 
changes in age and metallicity alter the shape of the residual 
spectrum in rather different ways; while changing age alters 
the overall slope of the residual continuum level, changing 
the metallicity introduces stronger steps into the residual 
contirmum level, especially at ~ 2850. This is perhaps not 
surprising since much of the detailed UV continuum shape 
produced by stellar populations around this age is domi- 
nated by extensive line blanketing. 
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Figure 5. residuals for single stellar population fits to the F4V star, IfR 4683. The models fitted have the following ages and 
metallicities: top: 3 Gyr, 0.01 Zq; middle; 3 Gyr, 0.2 Zq; bottom: 3 Gyr, Zq; left; 1 Gyr, 0.2 Zq; right; 7 Gyr, 0.2 Zq. Where 
{^Xdata < ^Xmodel)-) the quantity plotted is —{{ixdata " ^Xmodel) / '^)'^ ■ The central wavelengths of the mid-UV spectral features identified 
by Fanelli et al. (1992), lying within this wavelength range, are marked on the spectra with dashed lines. See §3.2 for discussion. 



For the sub-solar metallicity star, HR 4683, the spec- 
tral features identified by Fanelli et al. (1992), as marked on 
the plots, are remarkably well reproduced by the best-fitting 
model spectrum (figure 5, centre plot). The BL 2538 feature 
is least-well reproduced, although this is better fit by the 
best-fitting single stellar atmosphere model spectrum. One 
might expect that, as the spectral features of Fanelli et al. 
(1992) have been well studied, their flux would be well mod- 
elled. 

For the super-solar metallicity star, HR 4688, the Mgll 
2800 and BL 2538 features (offset in the model flux with 
respect to the central wavelengths of the features because of 
the binning resolution) are not well fitted by the best-fitting 
single stellar population model (figure 6, centre plot). It is 
unlikely that the poor modelling of these features arises from 
including post-MS contributions to the model flux, as they 
are not well modelled even by the best-fltting single stellar 
atmosphere model. The Mgll feature at 2800 A is better fit 
by the 0.5 model. This illustrates the dangers associated 
with attempts to determine age and metallicity using only 
a handful of absorption features. 



3.3 Mixed-metallicity models 

In figures 7 and 8, we show the results of fitting the mixed- 
metallicity model to the HST spectra of the F stars. Each 
figure shows the best-fitting model superimposed on the ob- 



served near-UV stellar spectrum, together with the resid- 
ual plots, with the central wavelengths of the mid-UV spec- 
tral features identified by Fanelli et al. (1992) marked by 
dashed lines. 

These are followed by contour plots of constant relative 
likelihood for the marginalised distribution of mean metal- 
licity with age. The contours contain 68.3%, 90% and 95.4% 
relative likelihood. Also included is a fourth panel showing 
how the fractional contribution (by baryonic mass) of the 
different metallicity components varies as a function of age. 
The results arc tabulated in table 3. 

This time, tiic interpolated results predict a mean 
metallicity for the F4V stax of 0.37 Z© at an age of 3 Gyr, 
while for the F8V star a mean metallicity of 1.52 Z©, at an 
age of 9 Gyr, is predicted. 

From the residual plot in figure, it can bo seen that 
the identified spectral features of HR 4683 are again well 
fitted by the model. For the super-solar metallicity stax, 
HR 4688, the absorption features BL 2538 and MGII 2800 
are much better reproduced by the mixed-metallicity model 
than by the singlo-motallicity models. 

Fitting this type of mixed-metallicity model arguably 
provides a fairer test of our ability to simultaneously deter- 
mine age and metallicity. Once again it can be seen that 
the derived metallicities and ages are well constrained, and 
the metallicities are in good agreement both with the com- 
pletely independent, optically-based results of Edvardsson et 
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Figure 6. residuals for single stellar population fits to the F8V star, HR 4688. The models fitted have the following ages and 
metallicities: top: 3 Gyr, 0.5 Zq; middle: 3 Gyr, 1.5 Zq; bottom: 3 Gyr, 2.5 Zq; left: 1 Gyr, 1.5 Zq; right: 7 Gyr, 1.5 Zq. Where 
{^Xdata < fxmociei)) ^he quantity plotted is —{(Jxdata " ^Xmodel) / '^)'^ ■ The central wavelengths of the mid-UV spectral features identified 
by Fanelli et al. (1992), lying within this wavelength range, are marked on the spectra with dashed lines. See §3.2 for discussion. 



al., and the results of the previous two subsections, although, 
in this case, the predicted age of HR 4688 is older than might 
be expected. However, it can be seen in the fourth plot of 
figure 8 that where a single metallicity dominates (at an age 
of ~ 3 Gyr) the results arc consistent with the results of the 
single-metallicity stellar population fits. 



4 RED GALAXIES AT Z ~ 1.5 
4.1 Single-metallicity models 

In figures 9 and 10, we show the results of fitting each of the 
seven single-metallicity models in turn to the Keck spectra 
of the red mjy radio galaxies LBDS 53W069 (Dunlop 1999, 
Dey et al. 2002, in preparation) and LBDS 53W091 (Dunlop 
et al. 1996; Spinrad ct al. 1997) respectively. The first plot 
in each figure shows the best-fitting model spectrum super- 
imposed on the rest-frame near-UV galaxy spectrum. The 
second plot shows contours of constant reduced {xi) a 
function of metallicity and age. The contour levels represent 
68%, 90% and 95.4% confidence limits. The results are pre- 
sented in table 2. An age of 3 Gyr, with a solax metallicity is 
derived for the galaxy LBDS 53W091 (z = 1.55). For LBDS 
53W069 (z = 1.43), ages greater than 6 Gyr are discarded, 
as they are not possible at z ~ 1.5 in currently-accepted 
cosmological models. The best fit age is then 3 Gyr, at a 
metallicity 1.5 Z0. 



In both LBDS 53W069 and 53W091, age and metaUicity 
have been simultaneously remarkably well constrained (see 
figures 9 and 10), with a predicted age for 53W091 between 
2.7 and 3.5 Gyr with 95.4% confidence, and for 53W069, be- 
tween 2.3 and 4.6 Gyr at the same level of confidence. Thus, 
under the assumption of a single metallicity for the observed 
stellar populations in these galaxies, it is clear that allowing 
metallicity to vary as a free parameter does not alter the 
conclusion of Dunlop (1999), that both these galaxies are 
~ 3 — 4 Gyr old, by any more than ~ 1 Gyr, even when, 
as in the case of LBDS 53W069, the predicted metallicity is 
super-solar. 

The results confirm the validity of previous arguments 
in favour of assuming roughly solar metallicity when mod- 
elling LBDS 53W091 (Spinrad ct al. 1997). 

We next investigate the robustness of these results by 
fitting the mixed-metallicity model to the galaxy spectra. 

4.2 Mixed-metallicity models 

In figures 10 and 11, we show the results of fitting the mixed- 
metallicity model to the Keck spectra of the LBDS galaxies. 
As in figures 7 and 8, each figure shows the best-fitting model 
superimposed on the rest-frame near-UV spectrum, followed 
by contour plots of constant as a function of mean metal- 
licity and age. The contours contain 68.3%, 90% and 95.4% 
relative likelihood. Again, a third panel is included show- 
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Figure 7. Top left: the best-fitting, seven-component, mixed-metallicity model (3 Gyr, 0.37 Zq, thin line - see §3.3 for details) to 
the spectrum of the F4V star, HR 4683 (thick line). Bottom left: residuals. Where {ixdata < fAmodeOi quantity plotted is 
~i{^\data ~ ^Xmodel)/'^)^ ■ The Central wavelengths of the mid-UV spectral features identified by Fanelli et al. (1992), lying within this 
wavelength range, are marked on the spectra with dashed lines. Top right: contour plots of constant relative likelihood, for the marginalised 
distribution of mean metallicity with age. The contours contain 68.3%, 90% and 95.4% relative likelihood. The shaded regions in the 
likelihood plots represent the estimates of Edvardsson et al., with estimated errors ± 30% in age and ± 0.10 dex in metallicity. Bottom 
right: fractional contributions (by baryonic mass) to the mixed metallicity model of the different metallicity components as a function of 
age , i.e. 0.01 Zq (thin solid), 0.2 Zq (thin dashed), 0.5 Zq (dot-dash), Zq (dash-dot-dot-dot), 1.5 Zq (dotted), 2.5 Zq (soHd) and 5.0 Zq 
(dashed) and 2.5 Zq (solid). See §3.3 for discussion, and table 3 for values of xi- 



object 



mean Z 
/Z© 



Best fit age 
/Gyr 



95.4% age Umits 
/Gyr 



F4V 
F8V 
53W069 
53W091 



0.37 
1.52 
3.78 
1.00 



2.7 - 3.4 
8.0 - 11.5 

> 4.5 
2.5 - 4.7 



1.47 
1.41 
0.78 
1.25 



Table 3. The results of fitting the near-ultraviolet spectra of the 
F-stars and the LBDS galaxies with the multi-component, mixed- 
metallicity model. Again, the fits are acceptable, although the fits 
to the stellar spectra are less good than the fits to the galaxies, as 
one might expect from fitting stellar population models to single 
stellar spectra. 



ing how the fractional contribution (by baryonic mass) of 
the different metallicity components varies as a function of 
age. The results are presented in table 3. Again, ages greater 
than 6 Gyr are rejected as impossible in currently accepted 
cosmological models. 



This time the best-fit mean metallicity of 53W069 is 
3.78 Zq. Only old ages (> 4.5 Gyr) are allowed, even at 
high metallicity. The third plot of figure 11 shows that the 
mixed-metallicity model is dominated by contributions from 
the high metallicity (2.5 and 5.0 Zq) models. This is a very 
different result from the single-metallicity fit, where the age 
is constrained at the 95.4 % level to between ^ 1.9 — 3.2 
Gyr, at less than half the metallicity here. 

Despite the relatively high value of best-fitting mass- 
weighted metallicity, the interesting thing is that the inferred 
age of LBDS 53W069 remains old. The second panel in figure 
11 demonstrates that, while the age is not constrained, ages 
< 4.5 Gyr are prohibited at the 95.4% confidence level. 

For 53W091, the mean metallicity derived from this 
mixed-metallicity modelling is reasonably well constrained, 
and is in fact precisely solar. Not surprisingly, therefore, the 
best-fit age remains ~ 3 Gyr (c.f. table 2), with ages < 2.5 
Gyr formally rejected at the 95.4% confidence level. Again, 
ages > 6 Gyr have been excluded on cosmological grounds. 

Is this lower age limit on 53W091 robust, or are better 
fits possible at lower ages with even higher metallicities than 
those investigated? For this galaxy, the bottom plot in fig- 
ure 12 shows that the highest metallicity component is only 
dominant at very young ages, where the statistical quality 
of the fit is completely unacceptable. A sufficient region of 
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Figure 8. Top left: the seven-component, mixed-metallicity model (9 Gyr, 1.52 Zq, thin line, see §3.3 for details) fitted to the spectrum 
of the F8V star, HR 4688 (thick line), which lies closest to the parameters corresponding to the minimum interpolated xi- Bottom 
left: residuals. Where {fxdata < ^Xmodel)^ the quantity plotted is —{{f\data ~ ^Xmodel)/'^)^ ■ The central wavelengths of the mid-UV 
spectral features identified by Fanelli et al. (1992), lying within this wavelength range, are marked on the spectra with dashed lines. Top 
right: contour plots of constant relative likelihood, for the marginalised distribution of mean metallicity with age. The contours contain 
68.3%, 90% and 95.4% relative likelihood. The shaded regions in the plots represent the estimates of Edvardsson et al., with estimated 
errors ± 30% in age and ± 0.10 dex in metallicity. Bottom right: fractional contributions (by baryonic mass) to the mixed metallicity 
model of the different metallicity components as a function of age , i.e. 0.01 Zq (thin solid), 0.2 Zq (thin dashed), 0.5 Zq (dot-dash), Zq 
(dash-dot-dot-dot), 1.5 Zq (dotted), 2.5 Zq (solid) and 5.0 Zq (dashed). See §3.3 for discussion, and table 3 for values of xi- 



parameter space has been explored for the age limit > 2.5 
Gyr to be regarded as robust. Again, it should be noted that 
where a single metallicity component dominates the fit, the 
results are consistent with the single-metallicity model fits. 

4.3 MOPED results 

In all of the above, we have assumed that the galaxies have 
formed the bulk of their stellar populations, of whatever 
metallicity, in a single coeval starburst. To check the extent 
to which this simplification is appropriate, it is clearly de- 
sirable to explore a larger parameter space, in which star 
formation history is allowed to vary freely. We have there- 
fore refitted the galaxy spectra using the MOPED algorithm 
(Heavens et al., 2000). In this approach, star formation is di- 
vided into several bins in time; the bins have equal width 
in log space, and allow widths small enough to account for 
the recent formation of stars in giant molecular clouds (10^ 
years). The height of each bin, allowed to fioat as a free 
parameter, represents the mass turned into stars at the cor- 
responding epoch. The metallicity of each bin is allowed to 
take any value from 0.01 to 5 Zq. In this way, we have the 
most model-independent prescription for describing a stel- 
lar population which is a composite of simple stellar pop- 
ulations with arbitrary metallicity, age and star formation 



rate. Since it is obvious that we cannot explore a parameter 
space with 1000 variables, we constrain ourselves to 12 bins 
in star formation history. With the help of MOPED (Heav- 
ens, Jimenez & Lahav, 2000, Reichardt, Jimenez & Heavens, 
2001), we are able to explore this huge (25 dimensional) pa- 
rameter space very efficiently, in only a few seconds. 

Figures 13 and 14 show the best fitting value found 
searching the hyper- likelihood surface*. The results for 
53W069 yield a population dominated by an old (4.25 — 7.72 
Gyr), high-metallicity (3.72 Zq) component, together with 
a small (7.0% by baryonic mass), younger component with 
a metallicity of 0.13 Zq. For 53W091, we can see an old 
component (2.35 — 4.25 Gyr) with a metallicity of 0.77 Zq 
and negligible contributions from other populations, in good 
agreement with our previous findings. This, therefore, con- 
firms that for both 53W091 and 53W069, the ages of the 
dominant populations remain old and apparently coeval, 
even when the star formation history is allowed to float 
freely. Unsurprisingly, therefore, the quality of the best fits 
achieved by MOPED are no better than those achieved with 

* The hyper-likelihood surface is so large in this case that not 
all of it can be explored at once. Instead we chose 2000 random 
starting points to find the global minimum, which is shown in 
figures 13 and 14, for 53W069 and 53W091 respectively. 
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Figure 9. Top: the best fit single metaUicity model (thin line, 
1.5 Zq model at 3 Gyr, after rejecting ages > 6 Gyr as cosmo- 
logically unreasonable) to the spectrum of the z = 1.43 galaxy, 
LBDS 53W069 (thick line). The central wavelengths of the mid- 
UV spectral features identified by Fanelli et al. (1992) are marked 
on the spectra with dashed lines. Bottom: contours of constant 
xt as a function of metallicity and age after interpolation. The 
contour levels represent 68%, 90% and 95.4% confidence. See §4.1 
for discussion, and table 2 for values of xi- 




age / Gyr 

Figure 10. Top: the best fit single metallicity model (thin line, 
Zq model at 3 Gyr) to the spectrum of the z = 1.55 galaxy, 
LBDS 53W091 (thick line). The central wavelengths of the mid- 
UV spectral features identified by Fanelli et al. (1992) are marked 
on the spectra with dashed lines. Bottom: contours of constant 
as a function of metallicity and age after interpolation. The 
contour levels represent 68%, 90% and 95.4% confidence. See §4.1 
for discussion, and table 2 for values of xi- 



the single burst models explored in the earlier section of this 
paper. 

It should be noted that the MOPED algorithm includes 
a parameter which allows for reddening of the spectra by 
dust. For both 53W091 and 53W069, the MOPED fits find 
a negligible contribution from dust. As the MOPED results 
are in good agreement with the results from the mixed- 
metallicity model spectral fits, we are justified in ignoring 
dust in the spectral fits. 



5 CONCLUSION 

The main conclusions of this work can be summarized as 
follows. 

First, we have demonstrated that, when fitted to high- 
quality, near-UV spectra of F stars, the synthetic near- 
UV spectra produced by the variable-metallicity models of 
Jimenez et al. (1998) are sufficiently accurate to extract well- 
constrained estimates of metallicity in accord with the val- 
ues determined from the study of well-established optical 
indices. This suggests that it should, in principle, be possi- 
ble to lift the age-metallicity degeneracy when fitting such 
models to high-quality optical spectra of high-redshift galax- 
ies. 

Second, fitting these variability-metallicity models to 



the Keck optical (rest-frame near-UV) spectrum of the red, 
z ~ 1.5 galaxy 53W091 we find that it is dominated by a 3- 
Gyr stellar population of solar metallicity, in good agreement 
with the original conclusions of Dunlop et al. (1996). 

Third, the fit to the even redder z ~ 1.5 galaxy 53W069 
has been able to constrain the metallicity well, but the age 
has not been convincingly determined. Given the success of 
the fits to the other spectra, where the old, 5 Z© model does 
not make any significant contribution, this is most likely to 
represent some uncertainty in the very highest metallicity 
model at old ages. However, better quality UV data and 
/ or the addition of rest-frame optical data may serve to 
resolve this problem. 

The main finding of importance for cosmology is that, 
even after this exploration of the effects of varying metallic- 
ity, the best-bet age of the oldest known galaxies at 2; ~ 1.5 
remains 3 Gyr, with ages younger than 2 Gyr now more 
strongly excluded than before. 

For Ho — 70kms~^Mpc~^ an Einstein-de Sitter uni- 
verse is younger than 3 Gyr at z ~ 1.5, and even the hard 
lower limit for coeval star formation of 2 Gyr translates into 
star-formation activity prior to z = 8. For a fiat universe 
with flm ~ 0.3 and JIa = 0.7, the best-bet age of 3 Gyr 
translates into z = 5, while the hard lower limit of > 2 
Gyr translates into 2: > 3. The added flexibility in the star 
formation history investigated using the MOPED algorithm 
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Figure 11. Top: the best-fitting, seven-component, mixed- 
metallicity model (6 Gyr, 3.78 Zq, thin Hne) to the spectrum 
of the high-redshift (z = 1.43) galaxy LBDS 53W069 (thick line). 
The central wavelengths of the mid-UV spectral features iden- 
tified by Fanelli et al. (1992) are marked on the spectra with 
dashed lines. Middle: contour plots of constant relative likeli- 
hood, for the marginalised distribution of mean metallicity with 
age. The contours contain 68.3%, 90% and 95.4% relative like- 
lihood. Bottom: fractional contributions (by baryonic mass) to 
the mixed-metallicity model of the different metallicity compo- 
nents as a function of age , i.e. 0.01 Zq (thin solid), 0.2 Zq (thin 
dashed), 0.5 Zq (dot-dash), Zq (dash-dot-dot-dot), 1.5 Zq (dot- 
ted), 2.5 Zg (solid) and 5.0 Zq (dashed). See §4.2 for discussion, 
and table 3 for values of xi- 




age / Gyr 

Figure 12. Top: the best-fitting, six-component, mixed- 
metallicity model (3 Gyr, Zq, thin line) to the spectrum of 
the high-redshift (z = 1.55) galaxy LBDS 53W091 (thick fine). 
The central wavelengths of the mid-UV spectral features iden- 
tified by Fanelli et al. (1992) are marked on the spectra with 
dashed lines. Middle: contour plots of constant relative likeli- 
hood, for the marginalised distribution of mean metallicity with 
age. The contours contain 68.3%, 90% and 95.4% relative like- 
lihood. Bottom: fractional contributions (by baryonic mass) to 
the mixed-metallicity model of the different metallicity compo- 
nents as a function of age, i.e. 0.01 Zq (thin solid), 0.2 Zq (thin 
dashed), 0.5 Zq (dot-dash), Zq (dash-dot-dot-dot), 1.5 Zq (dot- 
ted), 2.5 Zq (solid) and 5.0 Zq (dashed). See §4.2 for discussion, 
and table 3 for values of xi- 



allows some later star formation, but the bulk of the stars 
(> 90%) must still be formed at high redshift. 

Our results therefore continue to argue strongly against 
an Einstein-de Sitter universe, and favour a A-dominated 
universe in which star formation in at least LBDS 53W091 
is completed somewhere in the range 2 = 3 — 5. Such a con- 
clusion is in accord with recent results from sub-mm surveys 
which suggest that star-formation activity in radio galaxies, 



and potentially all massive ellipticals is largely completed at 
2 > 3 (Archibald et al. 2001; Dunlop 2001). 
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Figure 13. Top: best-fitting model (thick line) to the spectrum 
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al., 2000). Bottom: The predicted stellar population components. 
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with an age 4.25 - 7.72 Gyr. 
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